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Abstract The Euler equation for an inviscid, incompressible fluid in a three- 
dimensional domain M in R 3 implies that the vorticity is a frozen-in field. This can 
be used to construct a symplectic structure on R x M. The normalized vorticity and 
the suspended velocity fields are Hamiltonian with the function t and the Bernoulli 
function, respectively. The symplectic structure incorporates the helicity conservation 
law as an identity. The infinitesimal dilation for the symplectic two-form can be 
interpreted as a current vector for the helicity. The symplectic dilation implies the 
existence of contact hypersurfaces. In particular, these include contact structures on 
the space of streamlines and the Bernoulli surfaces. The symplectic structure on Rx M 
can be realized as symplectisations of these through the Euler equation. 



1 Introduction 

The traditional theories for the problems of fluid flow are statistical and proba- 
bilistic in nature [JJ , [5] • The geometric descriptions of the Eulerian equations of 
hydrodynamics in the framework of infinite dimensional Lie groups and Hamil- 
tonian (Lie-Poisson) formulations [3]- [2] has been attracted the attention of 
scientists working on different aspects of fluid motions. On the other hand, the 
need for some geometry, more specifically, the necessity of the symplectic ge- 
ometry for tackling some open problems of the Euler equation was argued. The 
extension of methods of symplectic geometry to viscous flows which is described, 
in Eulerian coordinates, by the Navier-Stokes equation, was also questioned [2]. 
The bottom line of these arguments may be the question of incorporating those 
geometric structures appropriate to and available from the Eulerian equations 
into the traditional statistical descriptions of fluid flows [5]-[TT].[I].[2]. 

Meanwhile, the relevance of geometric structures obtainable from the Eule- 
rian equations to the statistical description of turbulence was discussed in Ref. 

The Lagrangian description, that is, the description by the trajectories of 
the velocity field of steady fluid motion in two and four dimensions was shown 
to have the structure of a completely integrable Hamiltonian system [ID]. [IT]. 

In [12] , starting solely from the Euler equation of an ideal fluid we obtained 
an infinitesimal symmetry, a symplectic structure on R x M for the suspended 



velocity field and the conservation law for the helicity density. In [T3] and [2], 
we extended these intrinsic geometric structures to the study of kinematical 
symmetries and generalized helicity conservation laws of three dimensional in- 
compressible flows. In this work, we shall further exploit the Euler equation of 
ideal fluids to obtain geometric structures relevant to a qualitative study of the 
Lagrangian description of motion. 

We shall introduce a current vector field governing the dynamics of the he- 
licity density. This helicity current generates the dilation for the symplectic 
structure. With reference to the Navier-Stokes equation we shall demonstrate 
that the information content of the dynamical equations can be represented by 
the helicity current. We shall also discuss some algebraic consequences of the 
symplectic dilation. We shall then construct contact structures on the space of 
streamlines and on the Bernoulli surfaces. We shall show that the symplectic 
structure on the space of trajectories is symplectisations of these contact struc- 
tures. For fluid dynamical content of this work we shall refer to Refs. [2].[15] 
and the necessary mathematical background can be found in Refs. [H)]-[23]. 



2 An infinitesimal symmetry 

We shall begin with the Euler equation of ideal fluids 

— + v Vv = -Vp (1) 

for the divergence-free velocity field v tangent to the boundary of a connected 
region M C R 3 and the pressure function p. The identity v • Vv = ^V|v| 2 — 
v x (V x v) can be used to bring the Euler equation (TTJ) into the Bernoulli's 
form 

<9v 

— - vx(Vxv) = Va (2) 

where the function —a = p + v 2 /2 is the Bernoulli function [11], also called to be 
the total (or stagnation) pressure [24] . In terms of the divergence-free vorticity 
field w = V x v Eq. @ gives 

— - V x (v x w) = . (3) 

It follows from the identity 

V x (v x w) = (w ■ V)v - (v • V)w + (V • w)v - (V • v)w (4) 
together with V • v = V • w = that Eq. ([3]) is equivalent to 

— + [v,w]=0, [v, w] = (v • V)w - (w V)v (5) 

which means that w is an infinitesimal timc-dcpcndcnt symmetry of the velocity 
field v. That is, the time-dependent transformations generated by w on M 
leaves the trajectories of the velocity field invariant. 
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3 A symplectic structure 



In R 3 a vector field w corresponds to a unique two-form 

lo = w • (dx A dx) = rfjvf(v • dx) (6) 

where we let g?m = dx • V be the exterior derivative on R 3 . If we replace d^i by 
the exterior derivative d = dtdt+ c?m on R x M and use the Euler equation ^ 
to solve the time derivative v, f of the velocity field we obtain the two-form 

Q e = w ■ (dx A dx) - (v x w + Va) • dx A dt , (7) 

on R x M. By construction, fi e is closed on the space of solutions of the Euler 
equation. In other words, the three-form 

dtt e = (V • w) • (dx • dx A dx) + [w )t - V x (v x w)] • dx A dx A dt (8) 

vanishes for divergence-free vector field w satisfying the Euler equation in the 
vorticity form ([3|). Thus, we have the closed two- form Sl e on R x M as an 
extension 

u> = c?m(v • dx) — > tt e — d(v ■ dx) mod Eq. |(5J) 

of lo on M by the Euler equation and this is induced by the invariant differential 
operators 

dja - * d = dtdt + cJm • 

Under the condition that the Bernoulli function a is not the constant func- 
tion the skew-symmetric matrix of the two-form f2 e has a non-vanishing deter- 
minant. Equivalently, this can be expressed as 

ifi e A n e = -w ■ Va dx ■ dx A dx A dt ^ (9) 

whenever a ^ constant. It seems to be necessary to assume in addition that 
w/0. However, for a realistic fluid the validity of this condition was discussed 
in Ref. [25] . We thus obtained from the Euler equation a non-degenerate closed 
two-form Q, e on R x M, that is, a symplectic structure [IH]-[ST],[S2],[TB . The 
non-zero four-form in Eq. ^ is the symplectic or the Liouville volume element 
on R x M. The symplectic two- form is exact 

Q e = d6 e mod Eq.^, 6 e = v • dx + ip(t)dt (10) 

where tp(t) is an arbitrary function. 

4 Hamiltonian vector fields 

The non-degeneracy of tt e means that given a one-form (3 = (3 a dx a on R x M 
with the local coordinates (x a ) = (x° — t,x) the equation 

i(x)(n e ) = p, (n e ) ab x b = p a (ii) 
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has a unique solution for the vector field X — X a d a = X°d t + X • V and vise 
verse. Here, i(X)(-) denotes the interior product or the contraction with the 
vector field X [18], (^ e )ab are the components of the skew-symmetric matrix 
of the symplectic two-form Q e in the given coordinates and we employ the 
summation over repeated indices. 

For (3 being an exact one-form, namely, (3 = df for an arbitrary smooth 
function / on R x M the solution 

X f = ^^[w ■ V/(^ + v ■ V) - Jw V + (Va x V/) • V] (12) 

of Eq. (fTTj) is said to be the Hamiltonian vector field for O e with the Hamiltonian 
function /. From the skew-symmetry of (Q e ) a b we have the conservation law 

o = i(x f )i(x f )(n e ) = {n e ) ab x a x b = x a i a = x°f. t + x • v/ = o (13) 

for the Hamiltonian function. We observe that Xf reduces to — (w • Va) _1 w • V 
for the function f = t which obviously satisfies Eq. (fl~3|) . In other words, the 
normalized vorticity field is Hamiltonian 

i(X t )(Sl e ) = dt , X t = -(w • a) _1 w • V (14) 

with the Hamiltonian function t. 

A time-dependent conserved function for the velocity field can be found 
again from the Euler equation. We recall that an energy consequence of the 
Eulcr equation follows by taking dot product of its Bernoulli form with the 
velocity field. The result is known as the Bernoulli equation [2].|24j 

|(^ 2 )-v-Va = (15) 

which implies that if the pressure p does not depend explicitly on time, the 
Bernoulli function a is a time-dependent conserved function 

da 1 

^ + v-Va = 0, -a(i,x) = p(x) + -« 2 (i, x) (16) 

along the trajectories of the velocity field. Using this as the Hamiltonian func- 
tion we can write the suspended velocity field d t + v on R x M as a Hamiltonian 
vector field 

i(d t + v)(Sl e ) = da. (17) 

The symplectic structure VL e on R x M induces a Lie algebraic structure on the 
space of smooth functions on R x M with the Poisson bracket 

{f,g} e ^ J^ 1 )^ - x f^ = MX f ,X g ) (18) 

= ^[|w-V/-fwV, + (V/xV|,).Va] (19) 

where (£l~ 1 ) ab are the components of the inverse of the matrix of the symplectic 
two-form and d/dt — d/dt + v • V is the convective derivative. 
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5 Conservation laws 



We shall now discuss the differential invariants of the velocity field and obtain 
the conservation law for the helicity as an identity involving the invariant forms. 
To this end, we recall that a differential p— form £ is said to be a relative invariant 
for a vector field X if there exist a (p — 1)— form £ such that 

c x (0 = i(x)d£ + di(x)i; = d( (20) 

where Cx{-) is the Lie derivative. If C, = 0, £ is said to be an absolute invariant 
[2T],[TS]. It immediately follows from this definition, Eq. (TT7| and dil e — that 
the symplectic two-form is an absolute invariant of the suspended velocity field. 
In fact, this is true for any Hamiltonian vector field, so is for the normalized 
vorticity field X t . For the canonical one- form we compute the Lie derivative 

C dt+V (9 e ) =da + d{v 2 + ^(i)) = d(f/>(t) + X -v 2 - p) (21) 

using the Hamilton's equation (jTTJ) for the first term in the definition (|20p. Thus, 
8 e is a relative invariant. Moreover, the derivation property of the Lie derivative 
implies the relative invariance 

£ dt+v {6 e An e ) = £ dt+v {8 e )An e + 6 e A£ dt+v (n e ) =d((i/j{t) + ^v 2 -p)n e ) (22) 

of the three-form 

9 e A fl e = 2H dx ■ dx A dx + 

{{v 2 +-0)w - 2Hv + v x Vq) • dx Adx Adt (23) 

where the scalar component, namely, the coefficient of the term dx ■ dx A dx, is 
the helicity density 

We-vVxv=-vw. (24) 
2 2 v ; 

The absolute invariance of the four-form f2 e Att e follows from the same argument. 

Moreover, from the right hand side of Eq. © we conclude that the absolute 

invariance of Q e A fl e is a statement for the conservation of the Liouville density 

d 

— (w Va)+v- V(w- Va) =0 (25) 

along the trajectories of the velocity field. 

Since f2 e is closed and d9 e — Sl e these differential invariants of the velocity 
field satisfy the relation 

d(0 e A Q e ) - fl e A fl e = (26) 

identically. Using the expression © for the second term and computing the 
derivative of the last term in the three-form (1231) as 



d[(v x Va) • dx A dx A dt] — w ■ Va dx ■ dx A dx A dt (27) 



5 



we obtain the divergence expression 



¥ + V.(Hv--(-^ P )w).0 



(28) 



for the conservation law of the helicity density TL. Note that the conserved 
flux, which is the vector field under the spatial divergence, is independent of 
the function ip which we have introduced arbitrarily. We can take it to be 
zero. Then, it immediately follows from Eqs. (|2Tj) and (|22|) that for v 2 /2 — p = 
constant, 9 e and 9 e A Q e become absolutely invariant. Moreover, from Eq. (128|) 
we see that on the level surfaces of the function v 2 /2 — p the density TL turns 
into a conserved function of the velocity field. 

We thus obtained the conservation law for the helicity density TL = v • w/2 
associated with the vorticity field w which was shown to be an infinitesimal kinc- 
matical symmetry. This means that, although, the conservation law is given as a 
divergence expression the density is a conserved quantity of the Lagrangian mo- 
tion. The symplectic structure f2 e has a prominant role in relating the vorticity 
field w to the helicity invariant. 

We conclude that, apart from the assumption on the pressure function, the 
Euler equation itself intrinsically contains a geometric structure in which we 
can relate symmetries, invariant differential forms and conservation laws. In 
the next section we shall examine the situation for the Navier-Stokes equation. 
We shall then continue our discussion of geometric structures available from the 
Euler equation. 

6 Viscous flow 

The Navier-Stokes equation for a viscous incompressible fluid in a bounded 
domain M C R 3 is 



where v is the kinematic viscosity [2 . As in the case of the Euler equation, Eq. 
(j2"9")) can be brought into the form 




Vp + zA7 2 v 



(29) 



— V x (v x w) = z/V w 



(30) 



by which the two-form 



fi„ = w • (dx A dx) — (v x w — z/V x w) ■ dx A dt 



(31) 



can be shown to be closed. For non-degeneracy we compute 



-fij, A Q„ = 2isH w dx • dx A dx A dt ^ 



(32) 



where the density 




1 1 2 

— w • V x w = — w • V v 

2 2 



(33) 



G 



is known as the vortical helicity [2J. For viscous flows V 2 v ^ and, for a 
realistic fluid we have w/0. Hence TL W ^ 0. This makes Vl^ non-degenerate. 
Sl v is also exact with the same canonical one-form 9 e 

tt v = d9 e mod Eq. (J29J) (34) 

which shows that the symplectic structure on the space of trajectories manifests 
the dynamical properties of the fluid. 

For an arbitrary smooth function / on R x M the Hamiltonian vector field 
Xf defined by the symplectic two- form f2„ is 

X f = ^[-w • V/(| + v) + f W + „V/ x V x w) • V] (35) 

which gives (l/2v7i w )-w ■ V for / = t. The suspended velocity field is, however, 
not even locally Hamiltonian. These are vector fields obtained from Eq. 
for closed non-exact one-forms (3. The symplectic two-form is invariant under 
the flows of locally Hamiltonian vector fields because Cx (Qu) — di(X)(fl I/ ) = 
d(3 = where we used the identity Cx = i{X)od + doi(X) for the Lie derivative 
and dVL u = 0. We find, however, that the one-form 

i(d t + v)(Q u ) = ^V 2 v • (dx - vdt) (36) 

obtained by contraction with the suspended velocity field is not even closed due 
to diffusive term. That means the symplectic two-form is not invariant under the 
Lagrangian motion. A consequence of this non-invariance or, equivalently, of the 
viscous diffusion is the non-conservation of the helicity density. Nevertheless, 
using the symplectic form (f3"Tj) and 9 e in the identity ([26| we get 

^ + V • (Wv + vw x V 2 v - iv 2 w) = -1vH w (37) 

as the evolution equation for H. 

These are properties of the Lagrangian motion reflected by the dynamical 
(or Eulerian) equations. We shall continue our discussion of the geometry of 
the Lagrangian motion in the framework of the Euler equation. We shall occa- 
sionally present the results for the Navier-Stokes equation in order to see how 
the dynamics is encoded into geometric objects. 



7 A symplectic dilation 

For our discussion of the geometric structure of the Lagrangian motion of the 
Euler flow a particularly interesting solution of Eq. (fTTj) is obtained when we 
let the one-form /? be the canonical one- form 9 e . The vector field J satisfying 
the equation 

i(J)(n e )=9 e (38) 
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can be uniquely determined to be 

J = — —\2H{d t + v ■ V) + (<tp(t) - v 2 )w • V + v x Va • V] (39) 

for the Euler equation and, for the Navier-Stokes equation we find 

J » = 7ri-[^(dt + v) + (p- \v 2 )w - i/v x V x w V] . (40) 
lvri w I 

It follows from Eq. (|38|) and d£L e — that J fulfills the condition 

c,,(n e ) = di(j)(n e ) = de e = n e (4i) 

of being an infinitesimal symplectic dilation for f2 e [26j . As a consequence of 
Eq. (|4ip and the derivation property of the Lie derivative we see that J expands 
the Liouville volume in Eq. |9]). J is also said to be the Liouville vector field of 
VL e [H]. 

We shall now give an interpretation of the symplectic dilation J in connection 
with the helicity conservation. We first compute its divergence with respect to 
the symplectic volume [18] . [19] . [2T] . [23] 

(di«(i e J)(k a n e ) = Cj{-n e a n e ) = o e a o e (42) 

where we used Eq. (141) . By the same equation we see that the second equality 
is the identity in Eq. (|26|) resulting in the helicity conservation. Thus, the 
conservation law (128j) can be written as the equation 

div n J-2 = (43) 

involving the symplectic-divergence of the symplectic dilation. With this inter- 
pretation we intend to call J to be the current associated with the helicity. The 
dynamical content of the helicity current can be revealed from a comparison of 
the symplectic structures obtained from the Navier-Stokes and the Euler equa- 
tions. The canonical one-forms are the same. So, the dynamics is encoded into 
the symplectic two-forms. They define the current vectors by Eq. (|38[) for the 
same canonical one-form. With this definition, the dynamical properties of the 
fluid, such as viscosity, become implicit in the helicity current. Thus, we can 
think of the pair (# e , J) as a geometric representative of the dynamics of fluid 
motion on the space of trajectories. 



8 Algebraic consequences 

The helicity current is not a Hamiltonian vector field. However, it takes a 
Hamiltonian vector field into a Hamiltonian vector field by its action via Lie 
derivative. To see this, we compute 

i([j,x f ])(n e ) ^ £Mx f )(n e ))~t(x f )(£j(n e )) (44) 

= -d(J(f)-f) (45) 
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where Eq. flU) is an identity [53], [TS] and we used Eq. (gTJ). Thus, [J,X/] is 
Hamiltonian with the function J(f) — /. Replacing Xf with [J, X/] in Eq. (|44[) 
and using Eq. we get 

iQJjJ.X/]])^) = £ J ( i ([J,X / ])(f7 e ))-z([J,X / ])(£ J (fi e )) (46) 
= -d(J(J(f))-2J(f) + f) (47) 

which is also Hamiltonian. Thus, by repeated applications of the Lie derivative 
with respect to the symplectic dilation J one can generate an infinite hierarchy 
of Hamiltonian vector fields (Cj) k (Xf) anchored to Xf. 

In particular, we let / = t so that X t = — (w • Va) _1 w • V and consider the 
Hamiltonian vector fields {Cj) k [X t ) for k > 0. The identity [25].[Tg] 

£[J,x t ] = £.J ° £x t - C Xt o Cj (48) 
evaluated on the suspended velocity field dt + v gives 

£[j,x t ) (d t + v) = -C Xt ( [J, d t +v}) (49) 

where we used [dt + v, X t ] = 0. The vector field [J, dt + v] is Hamiltonian with 
the function J (a) — a — tp — v 2 . So, we can express the right hand side of Eq. 
(|49|) as a Hamiltonian vector field. To find the Hamiltonian function we shall 
use the Lie algebra isomorphism 

[X f ,X g ] = -X {f , g} (50) 

defined by the symplectic structure fi e [23].|21j. Taking / = t and g — ip — v 2 
we compute the Poisson bracket 

{ t ^-v 2 } = X t (iJ-v 2 ) = ^^- (51) 

which gives the Hamiltonian function up to a constant. It follows from Eq. (|49p 
that on the level surfaces defined by the constant values of the Hamiltonian 
function (f5Tj) of Cx t {[J, d t + v]) we have [[J, X t ],d t + v] = 0. Since w ^ 0, this 
means that w lies on the two-parameter family p + av 2 /2 = b, a,b = constants, 
of surfaces. In this case, the hierarchy of Hamiltonian vector fields anchored 
to X t becomes infinitesimal symmetries of the velocity field. In particular, this 
family of surfaces includes the Bernoulli surfaces for a = 1. In other words, on 
the level sets —a — b of the Hamiltonian function, d t + v admits the infinite 
hierarchy (Cj) k (X t ) of time-dependent infinitesimal symmetries. 

9 Contact structures 

A contact structure on a three dimensional manifold is a field of non-integrable, 
two-dimensional hyperplanes in its tangent spaces. Locally, this may be de- 
scribed as the kernel of a one-form to satisfying u A dio ^ everywhere. The 
contact form u> determines a unique vector field E by the conditions 

i{E){uj) = l, i{E){duj)=Q (52) 
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which is called the Reeb vector field. [2"0]-[2"2"]. 

A geometric consequence of the existence of a symplectic dilation is that 
the canonical one-form 9 e can be related to a contact structure on some three 
dimensional domain in R x M because from Eq. (|41[) we have 

e e a de e = e e a n e = i{j){n e ) a n e = i(J)(-n e a n e ) ^ o . (53) 

According to the definition of Ref. [28 a hypersurfacc in a symplectic man- 
ifold admits a contact one-form if and only if there exists a symplectic dilation 
which is defined on its neighborhood and is transversal to the hypersurface. 
It is also remarked that such hypersurfaces arise as level sets of Hamiltonian 
functions of Hamiltonian vector fields [26] -[28], [22] . In this section we shall de- 
scribe two examples of contact hypersurfaces in R x M . In the next section we 
shall consider their relation with the symplectic structure fi e and discuss some 
physical significance. 

9.1 Spatial hypersurfaces 

First, we consider spatial hypersurfaces M c = {(t, x) GRxM\t = c = 
constant} as level sets of the function t. Recall that t is the Hamiltonian function 
for the normalized vorticity field X t . Then, the transversality condition 

i(J)(i(X t )(Q e )) = i(J)(dt) = f (54) 

w • Va 

for M c holds for non-vanishing values of the hclicity density. In other words, 
for H 7^ the helicity current is not contained in the tangent spaces to M c . In 
this case, the contact form on M c is obtained as follows. Let 

i;M c ->Rx M : (t = c, x) i-> (t, x) (55) 

be the inclusion of time slices M c into space-time. A function onJixM gives 
a function on M c when composed by i. This operation can then be extended to 
differential forms. If a = f a (x)dx a is a one-form on R x M its pull-back i*a to 
M c by the inclusion map i is defined to be 

l*CJ = i*(f a (x)dx a ) = l*{f a (x))l*{dx a ) = {f a O l){x)dl( X a ) (56) 

where we used the commutativity of the operators d and i* [TH]. In particular, 
we compute 

oj = i*9 e = v(i = c,x) • dx (57) 
duj = di*9 e = i*d6 e = i*Q e = w(t = c,x) ■ dx A dx (58) 

and it follows that w is a contact form on M c 

lo A d M uj = i*9 e A i*fl e = i*(9 e A fi e ) = 2H dx ■ dx A dx ^ (59) 
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provided the helicity density is non-zero. Note that as a result of the transver- 
sality condition we have i(J)(ui) = 0. In fact, i(J)(Q e ) = follows from the very 
definition of the symplectic dilation J. 

We shall now show that the investigation of the normalized vorticity field 
X t as a Hamiltonian system on R x M is equivalent to the study of the Reeb 
vector field on the level surfaces M c of the Hamiltonian function t. Since t o i is 
the constant function on time slices M c we have 

i*[i(X t )(fl e )) = i*(dt) = di*t = d(toi) = Q (60) 

as the pull-back of the Hamilton's equations for X t to the spatial hypersurfaces. 
On the other hand, for any vector field X on R x M we have the identity 

i*(i(x)(n e )) = i(i*(x))(i*a e ) = i{u{x)){d M u) (6i) 

where i^X denotes the push-forward of X to M c . That means, i*X is the 
pull-back of X by and hence is a vector field on M c [18] . Since the one 
dimensional kernel of duu (in the tangent spaces of M c ) is spanned by the 
Reeb vector field, Eq. (|60j) for X t is possible only if the push-forward i*Xt is 
proportional to the Reeb vector field of the contact structure on M c . In fact, it 
is easy to check that the vector field 

E(c,x) = (^w)(f = c,x) = (^^X t )(t = c,x) (62) 

satisfies the criteria in Eq. (|52"[) . 

To this end, we want to remark that the contact structure on spatial hyper- 
surfaces M c must be distinguished from similar geometric constructions on M 
obtained from the Euler equation 

v(x) x w(x) = -Va(x) (63) 

for the steady flow of incompressible fluid. In this latter case there is also 
a contact structure on M provided the (time-independent) helicity density is 
non-zero. However, the difference between the two cases is not merely the time- 
dependence of fields. They imply qualitatively different descriptions of the flows. 
For example, from Eq. ([63]) of the steady flow we obtain 

v Va = , w • Va = , [v, w] = (64) 

which means that the fields v and w span the tangent spaces of the (two dimen- 
sional) Bernoulli surfaces a = constant and their flow lines commute on these 
surfaces [20] . On the other hand, the pull-back of the unsteady Euler equation 
to the hypersurfaces M c gives 

v )t (£ = c,x) — v(t = c, x) x w(t = c,x) = Va(t = c,x) . (65) 

The qualitative analysis of this equation implies quite different and complicated 
results for the surfaces a(t = c, x) = constant as well as the expression for the 
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push-forward of the suspended velocity field to the contact hypersurfaces M c . 
The behaviour of the flows lines on M c are also different. For example, if we 
take the curl of Eq. fBSjl 

w t t(t = c, x) + [v(< = c, x), w(t — c, x)] = (66) 

we see that the flows of the velocity and the vorticity fields do not necessarily 
commute on M c . 

9.2 Bernoulli surfaces 

Let B C R x M be the level sets of the Bernoulli function a which is the 
Hamiltonian function for the suspended velocity field. Let i : B — > R x M be 
the inclusion. The transversality condition reads 

i{J){da) = ip(t) -v 2 ^0 (67) 

where we used the fact that a is conserved under the Lagrangian motion. Since 
ip is arbitrary, the transversality of the Bernoulli surfaces to the hclicity current 
is the same as the non-vanishing of the kinetic energy of the fluid. On B the 
Eulcr equation becomes 

— vxw = 0, a = constant . (68) 

dt K ' 

Obviously, the pull-back of the symplectic two-form 

i*VL e = w • (dx A dx) - v x w • dx A dt (69) 

to Bernoulli surfaces is degenerate. Its one-dimensional kernel in the tangent 
spaces of B is the span of the vorticity field. Since d t + v is also tangent to B the 
two-dimensional tangent hypersurfaces (in the three dimensional tangent spaces 
of B) on which i*Q, e is non-degenerate can be defined to be the complement of 
span{dt + v, w ■ V} in the tangent spaces of R x M. i*fl e is also exact on B 

i*fl e = d(v • dx) mod Eq. JBSJ (70) 

by the pulled-back Euler equation. So, the contact structure on the Bernoulli 
surfaces is defined by the time-dependent one-form 

uj t — v(t,x) • dx (71) 

whose derivative is the two-form in Eq. (|69p . The non-integrability condition 

uj t A dw t = 2H dx ■ dx A dx + (v 2 w - 2Hv) ■ dx A dx A dt (72) 

of the tangent hyperplanes defined as above requires cither the hclicity density 
or the kinetic energy to be non-zero. Recall that the non- vanishing of the kinetic 
energy is also required by the transversality condition. For contact one-form ui t 
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on the Bernoulli surfaces we can find the Reeb vector field up to an arbitrary 
function 

E = m(t,x)(d t +v) + n(t,x)w-'V , mv 2 + 2nH=l. (73) 

This arbitrariness is a manifestation of the fact that contrary to the case of 
spatial hypersurfaces the inclusion of the Bernoulli hypersurfaces into R x M 
are defined only implicitly. 



10 Symplectisation 

If u> is a contact form on M with local coordinates x then the two-form 

x) = e t (d M ^(x) + dt A uj(x)) = d(e t u(x)) (74) 

defines a symplectic structure on R x M and is called to be the symplectisation 
of uj. Here, t is the coordinate on R. Note that has dt as its symplectic 
dilation and that the Reeb vector field E on M defined by to is a Hamiltonian 
vector field for the symplectisation with the Hamiltonian function e* [20]. [22]. 

We can think of the above symplectisation to be induced by the map u) i— > e u> 
in which the time variable is introduced artificially. The geometric fluid dynam- 
ics provides an unusual but nevertheless a natural example of symplectisation. 
Recall that we obtain the contact one-form on time slices by pulling the canon- 
ical one-form 9 e back to M c by the inclusion. Conversely, the symplectisation 
of the contact structure on time slices follows from the inclusion map 

i : M c -> R X M : (t = c, x) i-> (i, x) . (75) 

In this case, the time variable t is introduced naturally by the action of the 
invariant differential operators. 

From a physical point of view the symplectisation of the time slices M c corre- 
sponds to the construction of trajectories of the velocity field from streamlines. 
These are solutions of the non-autonomous and autonomous equations 

^=v(*,x), ^=v(t = c,x(r)), (76) 

respectively. The solutions to the first equation onExM can be constructed 
by solving the autonomous system on M at each time t and then joining them 
by the inclusion 

i : (streamlines) >—> (trajectories) . (77) 

The symplectisation to R x M of the contact structures on M c means that the 
inclusion in Eq. (|77p for solutions of the differential equations extends to 

contact structure on \ / symplectic structure on 

the space of streamlines J \ the space of trajectories 

which covers the corresponding geometric structures. 
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The symplectisation of the contact structure on the Bernoulli surfaces may 
be given a similar interpretation in the language of the solutions of differential 
equations. In this case, it will be appropriate to consider the solutions of the 
Euler equation. The inclusion 



of the solutions of the Euler equation from the solutions of a homogeneous equa- 
tion on each Bernoulli surface. This may be interpreted as a nonlinear analog of 
the result in the theory of linear differential equations that the general solution 
of an inhomogeneous equation is the sum of the solution of its homogeneous 
part and the particular solution. 

11 Conclusions 

We obtained geometric and algebraic structures for the Lagrangian motion of 
an incompressible fluid implied by the Euler equation. The point of the paper 
is that the dynamical Eulerian equations of ideal fluid determine the geometry 
of the space of trajectories and can be used for qualitative investigations of the 
Lagrangian motion. 

The most intriguing result in this direction may be the introduction of the 
symplectic dilation. We showed that the dynamical equations can be represented 
geometrically by a pair consisting of the canonical one-form and the symplectic 
dilation. As an immediate consequence of the existence of a symplectic dila- 
tion we the evolution equation for the helicity density and studied the contact 
structures on the spatial hypersurfaces and on the Bernoulli surfaces. 

Since, the transformation generated by the symplectic dilation expands the 
Liouville volume we can argue that it generates the scaling transformation on the 
space of Eulerian fields. In spite of the general belief that a precise connection 
exists between the scaling transformations and helicity conservation law, we 
have not been able to find any source in the literature. 

We obtained the helicity conservation law as an identity among invariant 
differential forms associated with the vorticity field. This manifests a relation- 
ship between symmetries and conservation laws which is different from the one 
familiar, for example, from the Noether theorems. Being realized as a symmetry 
of the flow of an ideal fluid the vorticity flow turns out to be the underlying 
symmetry of the Lagrangian motion connected with the helicity conservation. 




( space — time ) 



implies the construction 
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